ALTERNATIVE JACOBI POLYNOMIALS 
AND ORTHOGONAL EXPONENTIALS 

VLADIMIR S. CHELYSHKOV* 

Abstract. Sequences of orthogonal polynomials that are alternative to the Jacobi polynomials on the interval 
[0, 1] are defined and their properties are established. An (a, /3)-parameterized system of orthogonal polynomials 
of the exponential function on the semi-axis [0, oo) is presented. Two subsystems of the alternative Jacobi polyno- 
mials, as well as orthogonal exponential polynomials are described. Two parameterized systems of discretely almost 
orthogonal functions on the interval [0, 1] are introduced. 
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1. Introduction. The concept of inverse orthogonalization of a sequence of functions 
emerged from the study of spectral methods for numerical simulation of near-wall turbu- 
lent flows |6]|3). Turbulence in the boundary layer near a fiat plate is an intensive physical 
phenomenon that exhibits an exponential decay of disturbances far from the wall. Thus, 
the sequence of exponential functions {e _fcj, }jj =1 , y <E [0, oo), being orthogonalized, can be 
employed for numerical simulation of the phenomenon. It was noticed that applying the al- 
ternative algorithm of inverse orthogonalization to the sequence of exponential functions is as 
constructive as applying the direct orthogonalization algorithm to the sequence of monomials 
{x k }^ =Q , x € [0, 1], and the bidirectional algorithm of orthogonalization was introduced in 
for defining systems of orthogonal exponential polynomials on the semi-axis. 

Later, the alternative Legendre polynomials were presented, and two associated quadra- 
ture rules were obtained . 

The alternative Jacobi polynomials are described in this paper. By taking advantage of 
the bidirectional orthogonalization we establish general properties of the polynomials. We 
show that the shifted Jacobi polynomials corresponding to the weight function with a trans- 
lated parameter are a special subset of the alternative Jacobi polynomials. Also, we consider 
an exceptional case of singular orthogonality; this feature takes place in some of the alterna- 
tive systems. 

We find that the alternative Jacobi polynomials in an exponential form are the system 
of functions on the semi-axis [0, oo) that has the same algorithmic capacity as the Jacobi 
polynomials have on the interval [—1, 1]. Such a deduction is illustrated with two special 
subsystems of the exponential polynomials possessing nice properties. 

The orthogonal exponential polynomials on the semi-axis may serve a useful purpose in 
approximation of a continuous function on a closed interval. Following [4], we introduce two 
parameterized systems of discretely almost orthogonal functions on the interval [0, 1]. The 
systems are constituted of the orthogonal exponential polynomials and unity. For the second 
system, supplementary computations are required to meet the definition of the functions, and 
properties of such a special construction have not been studied. 

To obtain these results we followed the classical theory of orthogonal polynomials. 
Different aspects of the theory that we were inspired by are represented in monographs 
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2. Construction of an (a, f3) -Parameterized Alternative Orthogonal Polynomials. 

Let n be a fixed whole number and 



vt fj \x) = {v^\ X )r k=n 

is the system of polynomials defined by Rodrigues' type formula as 



(2.1) 



- a -k-l( 1 _ x yf3 d n-k 



[n — k)\ dx T 
with a > — 1, /3 > — 1. From ( |2.2| it immediately follows that 



J- 



notation (a) m represents the falling factorial. The last formula contains the term of lowest 
order in x such that 



T(a + n + k + 2) k 
(n- k)\r{a + 2k + 2) X 



where T(z) is the gamma function. 
Lemma 2.1. For < k < I < n 



x a (l - xfv^ (x)x l dx = 0. 



(2.3) 



(2.4) 



Proof. Substituting (2.2 1 to (2.4 1 and integrating by parts we confirm (2.4 1. □ 
THEOREM 2.2. The alternative Jacobi polynomials V^l (x) sa ^ s fy 



x a {1 _ x f V ^\ X )V^\x)dx = h^S kl , (2.5) 



l nk 



1 T(a + n + k + 2)T{(3 + n- k + 1) 
+ 2k + l (n- k)\r(a + f3 + n + k + 2) 



(2.6) 



Here Ski is the Kronecker delta. 

Proof. Since orth ogon ality i s ve rifie d by Lemma 1.1, we only need to consider the case 
/ = k. Substituting (2.2 1 and (2.3 1 to (2.5 1 and integrating by parts we calculate 



we confirm the restrictions on a and (3 in (2.1 



Substituting k = and k = n to h , 
□ 

COROLLARY 2.3. For p £ N orthogonality relations \2.4\ - \2.6\ hold the invariance 



■(a-2p,/3) 



L nk ' 



and 



x p V. 



nk 



(x) 



,(a-2p,/3) 



(*). 



ALTERNATIVE JACOBI POLYNOMIALS 



3 



Polynomials Vp"£ (x) — x v V^ k (x) are orthogonal on the interval [0,1] with the 
weight function x a ~ 2p (l — x) 13 . 

Following the proof of Theorem 2.2 we also evaluate 



x f V ^\x)dx = r(q + fc + i)r03 + n-fc + i) 



k\ 



(n-k)\ T{a + /3 + n + 2)' 



Making use of \2.2) , Cauchy's integral formula for analytic functions, and reciprocal 
substitutions one can obtain the integral representation 

, „, 1 ( rr -l\a+P+n+2 r -(Q+/3+n+fe+2)/i _ \/3+n-k 

p^'fifx) = _ \- > / f \t Z A dz (2 7) 

1 nk W { 1 „-X\B I („ „.-l\ra-fc+l ' y ' 



2m (l-x- 1 )! 3 



(z-x- 1 ) 



where C is a closed contour encircling the point z — x 1 . 

2.1. Reciprocity. Being constructed by the procedure of inverse orthogonalization, poly- 
nomials V^' (x) can be represented in terms of the Jacobi polynomials. Below we derive 
two relationships of such kind. 

First we determine the relationship between "P^ '^(x) and polynomials from the se- 
quence 

P^ ) (x) = {P^'«(x)}^ = „, 
that are defined by the procedure of direct orthogonalization as 

sign (P^'«(l))= (-!)*-», (2.8) 



x«(l - xfP^\x)P^\x)dx = d^S kl , 



(2.9) 



1 



T(a + k + n + l)T{/3 + k-n + l) 



nk a + /3 + 2fc + l (k-n)T(a + fi + k + n + 1) 



(2.10) 



Since n is fixed, by verifying (2.8 I, (2.9 1, and (2.10 1, the polynomials P^' (x) can be 



immediately connected to a fixed set of the Jacobi polynomials P^f ' (x) Ifl2l . Precisely. 



P^ ) (x)=x»Pi^ 2 "'«(l-2x). 



(2.11) 



This formula can be used directly for describing properties of P^J (x), and one of the 
results that follow from ( |2.1 \\ is the integral representation 



27ri(l-x)/ 3 7 (z-x) fe -™+ 1 



dz. 



(2.12) 



Here Ci is a closed contour encircling the point z = x. Comparison of ( 2.7 1 and ( 2. 12 1 leads 
to 



(2.13) 
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Thus, relation of reciprocity (2.13 1 associates two sequences of polynomials that contain 
different powers x k - i.e., with < k < n and with n < k < 2n respectively. 

Let us now consider the relationship between 'P„ (x) and polynomials from the se- 
quence Pq ' (x) such that 



P^\x)=P^\l-2x). 



(2.14) 



Shifted Jacobi polynomials (2.14i are orthogonal on the interval [0, 1] with the weight func- 
tion x a (l — x)P . The polynomials can also be defined by the Rodrigues formula as 



Comparing ( 2.2 1 and ( 2. 15 1 we find that 



(2.15) 



(2.16) 



2.2. Recurrence Relations and Properties Involving Differentiation. Relation of reci- 
procity (2.13 1 and formula ( |2.1 \) , as well as (2.16 1 and ( 2.14| >, enable employing the classical 



results for establishing properties of V^l (x). 

Following ( |2.13| l and ( |2.11| i we obtain the three-term recurrence relation 

V<&»(x)=x n , 



VnX-M = + ^n)x n - x - (a + (3 + 2n + l)x n 



(n-k + l)(a + n + k + l)(a + 2k + 2)V^f\{x) = {a + 2k + 1) 



x [(a + 2k)(a + 2k + 2)x~ 1 - (a + 2n + 2){a + f3 + 2k + 1) - 2(n - k)(n - k + 1)] 

xV ( n l' P \x) - (a + + n + k + 2)08 + n - k){a + 2k)T ( n a k %(x) 
and the differential-difference relations 



(a + 2k + 2)x(l - x)^-V^f\x) = [ka + 2k(k + 1) - (to + n 2 + k l + 2n) 



dx 



xxr { n a /\x) -(a + p + n + k + 2)(/3 + n - k)xV^ p '{x), 



(«,/3) 



{a + 2k)x{l - x)—V { "f\x) = [-{a + 2k){a + k + 1) - f3{a + 2k) 



+((n + l)(a + /3 + n + l) + (a + k)(a + + k) + a + 2k)x]xV ( n l' 0) (x) 



+(n -k + l){a + n + k + l)xP^\(x) . 
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Following (2.16 1 and (2.14 1 we find the differentiation formula 



gv^\x) - kx-*r%»{z) =-(a + P + n + k + 2)V^f +1 \x) 
for k < n, and the differential equation 

x 2 (l - x)y"{x) + x[a + 2 - (a + /3 + 3)x]y'(x) 



(2.17) 



[fc(a + k + 1) - rc(a + (3 + n + 2)x]y(x) = 



that has a solution y = V^L (x). 

Let us consider the case k = 0. It follows from ( |2.16| l that 



a>-2, /?>-!. 



Thus, the shifted Jacobi polynomials with translated a, Pq^ +1 (x), 



(2.18) 



(2.19) 



are a subset of the 



alternative Jacobi polynomials vj^t {x). Relation (2.19 1 can also be directly traced from 



equation (2.18 1. For k — the degrees of the polynomial coefficients in (2.18 i are reduced 



by one, and this makes (x) the sequence with classical properties. 



of a ~P^fk^i(x) with an (a', /3'), and so dV^Z'^ (x) / dx do not belong to the original class 

of polynomials. Still, other properties of V^.'^ (x) described above are similar to properties 
of the classical orthogonal polynomials. 

3. Singular Orthogonality and Marginal Se que n ces. The restriction a > —1 for the 
the weight function w(x; a, j3) 
on a in (2.19 1. Let a' — {a : 



For k > polynomials (x) are less perfect. By considering property ( 2. 17 1, or by 

differentiating 



2. 18), we find that the first derivative of a V^.'^ (x) is not a constant multiple 



in (2.4 1, (2.5 i is stronger than the restriction 



x Q (l - 

—2 < a < —1}. In this exceptional case the polynomial 
T-^o ' p> (x) is not integrable with respect to the weight function w(x;a',(3) and therefore non- 



normalizable. However, it is consistent with orthogonality relations ( |2.4| i, ( |2.5| l for I > 0. 
Thus, "P^q '^(x) is a singular term of the sequence '^(x), which we call a marginal 
system. 

Assuming k > one can employ the sequences V^. (x) for approximation of a con- 
tinuous function /(x) on the interval [0, 1]. Such an approximation is effective near x = 
only if /(0) = 0. 

Below we describe two noteworthy orthogonal sequences that hold singular terms. 



3.1. A-Kind Orthogonal Polynomials. Considering (2.2 1 we find that P, ( l0 1,0) (x) are 



the shifted Legendre polynomials. This classical set of polynomials is a singular subset of the 
system 



A n (x) = T ( - lfi) (x), A n (x) = {A nk (x)}° k=n , neN. 



(3.1) 



The system A n (x) was introduced in [2|, and it is the first example of an inversely orthogo- 
nalized sequence of monomials {x k } k=0 . 

The polynomials A nk (x) obey the orthogonality relations 

fc = 0,l,...,n, i = l,2,...,n 



-A nk (x)A n i(x)dx = 



k + V 
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~ A^nk [x)dx — — , k — 1, ...,n. 
x k 



They can be calculated using the expansion 




FIG. 3.1. A-kind orthogonal polynomials: n = 5,k = 1 — 5. 

yU(*) = £(-iF ( n - j k )( n+ n k _ + k j )* k+ ^ * = 0,1 n, 

or following the three-term recurrence relation 

A nn (x) = x n , A^-iO) = (2n - l)x n ~ l - 2nx n , 

(2k + l)(n + k)(n -k + l)A„, fc _i(x) 

= 2fc[(2fc - l)(2fc + 1> _1 - 2(n 2 + k 2 + n)]A nk (x) 

-(2k -l)(n- k)(n + k + l).A„ )fe+1 (a;). 
They satisfy the differential-difference relations 



(2k+l)x(l-x)—A n . k (x) 
ax 



and 



= [2k 2 + 1 - (n 2 + k 2 + 2n)x]xA nk (x) - (n - k)(n + k + l)xA 7hk +i(x), 



(2k-l)x(l-x)—A nM (x) 
ax 



= [-k(2k - 1) + (n 2 + k 2 + n)x]xA nk (x) + (n + k)(n -k+ ^xA^k^x). 
We also find that y(x) = A n , k (x) is a solutions to the differential equation 

x 2 (l — x)y" (x) + x(l — 2x)y (x) — [k 2 — n(n + l)x]y(x) = 0. 
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3.2. T-Kind Orthogonal Polynomials. From (2.2 1 it follows that n\/(n - l/2)„ x 

V n0 2 ' 2 (x) are the (shifted) Chebyshev polynomials. This classical set of highly demanded 
in approximation theory orthogonal polynomials is a singular subset of the system of polyno- 
mials 



T„(x) = {T nfc (x)}° = „, 7 nk (x) = c nk V { nk ^ 5) (x), 

where c nk = (n — fe)!/(n — k — l/2) n _ k and neN. 
The polynomials obey the orthogonality relation 



(3.2) 



1 T nk (x)7 n i(x) ^ = 1 (2n-fc-Q!! (2n + fc + /-l)!! 

o xy/x(l - x) X 2{k + l)-l{2n-k-l-l)\\(2n + k + l-2)\\ 



irS k i, (3.3) 



and 



fc = 0,l,...,n, i = l,2,...,n. 



1 T nfe (x) (2fc~3)!! 

(2fc)!! ' fe = 1 ' 2 - n - 



/0 Xa/x(1 - l) 

They can be calculated using the three-term recurrence relation 




FIG. 3.2. T-kind orthogonal polynomials: n = 5, k = 1 — 5. 



T nn (x) = x", T rMi _i(x) - (4n - 3)x"~ 1 - (4n - 2)x n , 



[2(n + fc) - l][2(n - ft) + l](4fc + l)T n , fc _i(a;) 



(4ft - l)[(4fc - 3)(4fc + ^x" 1 - 2(4n 2 + 4fc 2 - 2k - l)]T nfc (x) 



-4(n -*)(»+ fc)(4ft - 3)T„, fe+1 (x). 
We also find that y(x) = 7 n<k (x) is a solution to the differential equation 

x 2 (l - x)y"{x) + x(l/2 - x)?/(x) - [k(k - 1/2) - n 2 x]y(x) = 0. 
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4. An (a, (3) -Parameterized Orthogonal Exponential Polynomials on the Semi- Axis. 

The classical orthogonal polynomials have a property of discrete orthogonality and the alter- 
native orthogonal polynomials have it also (see, for example, Let us consider the 
sequence 

exp(-fct), k e N. (4.1) 



Since 1 is not the term of the sequence, direct orthogonalization of (4. 1 1 on the semi-axis with 
a chosen weight function does not result in Gauss' type quadrature and, consequently, in the 
discrete orthogonality of the sequence. Thus, inverse orthogonalization is the only procedure 
for constructing the exponential polynomials possessing such a property (see, for example, 
0). This makes the set of orthogonal functions introduced below a unique system. 

Let £^'P\t) = (i)}fc_„ be the polynomials of exponential function defined as 

4l'* ) (t)=-Pir m (e- t )- (4-2) 



From (4.2 1 it immediately follows that S^'^' (t) is the orthogonal system, such that 



e~ at < 



1 T(a + n + k + l)T((3 + n-k + l) 

-Oki (4.3) 



a + 2k (n-k)W(a + P + n + k + l) 



for a > —1, P > -1. 

Properties of Sjft'^ (t) that are of interest for computations can be easily derived from the 



results of previous sections. In particular, ( 4.3 1 leads to the three-term recurrence relationship, 



which sets the associated orthogonal function £^^\t). Since £^'^\t) is not integrable 
on the semi-axis with the weight function 1, we do not include it in '^(t). Functions 
{gi^'P) (j S^'^^J (t) form the system with nice properties, because the zeros of £^^\t) 
are the abscissas of Gauss' type quadratures for the exponential functions on the semi-axis. 

In the next subsection we present two systems that are the exponential form of polyno- 
mials A nk (x) and 7 nk (x). 



4.1. Two Systems. Let a = p = 0. Following \A.2\ and ( |3.1| l we define 

£ nk {t) =A nk {e~ t ). 

The orthogonal exponential polynomials £ n k{t) were first introduced in [2]. Since they are 
of interest for approximation of a function in L 2 [0, oo), we repeat some properties of A nk (x) 
in the exponential form as follows 

1 

£nk{t)£ n i(t)dt = —Ski, k,l = l,...,n, 
2k 







J £ nk (t)dt = -. 



k = 1, n, 



'Formula (3.5) on page 21 in [5] have to be finalized as 

2 1 

w s 



n(n + 2) x1V nl (x s )V' n0 (x a ) 
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£nn(t) = e~ nt , £ n , n -i(t) = (2n - l)e- ( "- 1)4 - 2ne~ n \ 
(2k + l)(n + fc)(n - k + l)£ n , fc _i(t) 
= 2fc[(2fc - l)(2fc + l)e* - 2(n 2 + fc 2 + n)]£ nfe (t) 
-(2fc - l)(n - fc)(n + + l)£^ )Jk+ i(t), fc = n - 1, 1, 

— [f„, fc -i(t) +^, fc (t)] = -(fc - l)£ n , k -i(t) + k£ nk {t). 



Let a = (3 = -1/2. Following (4.2 1 and ( 3.2 1 we define 



For convenience, we also repeat some properties of T„fe (a;) in the exponential form as follows 



°° .u_ 1 (2(n -*))!! (2(n + fc)-l)H 



= di = — 1 ~ 7TT7 ~TZT~, ~ — 7T<5fci, fc, / = !,..., 71, 



o ^-'(l-e-*) 4fc-l(2(n-fc)-l)!!(2(n + fc)-2)!! 



l" , £lk{t) ^=^^=#2^, fc = l, 

= e-"*, = (4n - 3)e-(- 1 ) t - (4n - 2)e~"*, 

[2(n + fc) - l][2(n - fc) + l](4fc + l)^ fc _ x (*) 
= (4fc - l)[(4fc - 3)(4fc + l)e* - 2(4n 2 + 4fc 2 - 2fc - l)]£&(t) 

-4(n - fe)(n + fc)(4fc - 3)££ fe+1 (t), fc = n - 1, 1. 

The zeros of i s eas y t0 calculate, and this advantage can be a motivation for approxi- 

mation of a function on the semi-axis by £% k (t). 

4.2. Discretely Almost Orthogonal Functions on the Interval [0, 1]. Employing ex- 
ponential polynomials £^1'^ (t) can be of interest for approximation of a function on the 
interval [0, 1]. 

One may try to apply the system 

£ { n a >"\t) = {l,£^\tm =1 , £^\t) = <?&«(A&0 • t). 

as the basis functions. Here, the scale parameter aI" ' 3 ' is the maximum non-trivial zero of 
the zeros \ { " k ' 0) of the function £^\t) 0. 

Another way is the use of parameters (a, (3). Below we describe a formal construction 
of functions that involves computation of (a n ,/3 n ) for desirable distribution of the zeros 

A lfe"' /3 " ) on the interval (0, 1]. 
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First, by considering j3 = uia, we introduce parameter cu that sets the abscissa of the 
weight function maximum, t = ln(l + lu). Then, by subjecting choice of a n to the require- 
ments 



7r 



<1 and 7 n = maxA| 1 Q „"'"" a " ) , (4.4) 



we introduce parameter 7„. In general, the equality in (4.4 1 can be reached if a n ,/3 n G R; the 
inequality holds if a n , f3 n are selected from the set of whole numbers or from a set of rational 
numbers of interest. Accordingly, for j n = 1 we define the system of functions 

and, for -f n < 1, 

Similarly, we define the sequences Z^ n \t) and Z^"\t). Choice of u n sets distributions of 
the zeros of the functions under consideration, and Z^ n \l) = Z^ n \l) = 0. The system 

of functions (^Z^ n U Z^ n ^ (t) is easier to construct, in particular for a n , (3 n G N. 

An interpolating sequence of functions corresponding to approximation by Z-functions 
with n = 2, u-2 = 1 and ct2 , /?2 G N was employed in H for solving a stiff initial value 
problem. 

5. Conclusion. Making use of the algorithm of bidirectional orthogonalization we ob- 
tained two systems of orthogonal functions, i.e. the alternative Jacobi polynomials V^L [x) 

and the orthogonal exponential polynomials S^'^Xt). Algorithmic properties of the systems 
are very similar to the properties of the Jacobi polynomials. This indicates that each of the 
systems may be a convenient tool for applied analysis. 

In addition, we described two marginal sequences A n k(x) and 7 n k(x), as well as their 
exponential counterparts £ n k(t) and £^ k (t). Similarly, other important special systems of 
polynomials can be derived from the alternative Jacobi polynomials. 

Adapted systems of functions £ n ' (t) and Z^"\t) and their extensions that are ap- 
propriate for differentiation may be of interest for /zp-approximation on a finite interval. 

It should be mentioned that there are regular orthogonal exponential polynomials that are 
often used for approximation on the semi-axis, but usually they are not considered as special 
functions sui generis (say, in ifTTI ). Description of these orthogonal exponentials can be found 
in (8). Also, a special orthogonal sequence of exponentials that are defined on the semi-axis 
in a particularly appealing form is constructed in flj. 

The concept that is developed in this paper is quite consistent with the maxim of C.G.J. Ja- 
cobi: " Invert, always invert ". Many years ago we unknowingly subscribed to this dictum, 
and now we recognize its notability 1 9 1 . 
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